Is doped BaBiOs a conventional superconductor? 
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We report density functional calculations based on local density approximation (LDA) of the 
properties of doped barium bismuthates. Using linear-response approach developed in the frame- 
work of the linear mufiin-tin-orbital method the phonon spectrum of the Bao.6Ko.4Bi03 system is 
calculated and is compared with the results of the neutron diffraction measurements. The effect of 
doping in the calculation is modelled by the virtual crystal and mass approximations. The electron- 
phonon coupling constant A is then evaluated for a grid of phonon wave-vectors using the change 
in the potential due to phonon distortion found self-consistently. A large coupling of the electrons 
to the bond-stretching oxygen vibrations and especially to the breathing-like vibrations is estab- 
lished. Also, a strongly anharmonic potential well is found for tilting-like motions of the oxygen 
octahedra. This mode is not coupled to the electrons to linear order in the displacements, therefore 
an anharmonic contribution to A is estimated using frozen-phonon method. Our total (harmonic 
plus anharmonic) A is found to be 0.34. This is too small to explain high-temperature supercon- 
ductivity in Bao.6Ko.4Bi03 within the standard mechanism. Finally, based on standard LDA and 
LDA-I-U like calculations, a number of properties of pure BaBiOs such as tilting of the octahedra, 
breathing distortion, charge disproportionation and semiconducting energy gap value is evaluated 
and discussed in connection with the negative U extended Hubbard model frequently applied to this 
compound. 



I. INTRODUCTION. 

Since the discovery of s«eerconductivity at Tc ^30K 
in Bai_a;Ka;Bi03 (BKB0)EJI3, and from earher studies of 
BaPbi_:j;Bia;03 (BPBO) system with ~13K, there is a 
fundamental question whether the conventional phonon 
mediated pairing mechanism is operative in these high- 
Tc superconductors (HTSC). Doped barium bismuthates 
are different from the HTSC cuprateaj, since no anti- 
ferromagnetic ordering exists for the parent compound 
BaBiOs. This seriously doubts that strong electron cor- 
relations exist and are responsible for the pairing. Simple 
cubic superconducting phase makes BKBO and BPBO 
similar to the isotropic low Tc superconductors. How- 
ever, there is a number of features which makes doped 
bismuthates similar to the cuprates. Both systems are 
perovskite oxide superconductors with surprisingly low 
density of states at the Fermi level. This can hardly give 
high transition temperatures for the BCS-like supercon- 
ductors. As high-Tc cuprates are originated from antifer- 
romagnetic insulators, the parent BaBiOs compound is 
a charge-density wave (CDW) insulator in which oxy- 
gen octahedra around the Bi ions exhibit alternating 
breathing~in and breathing-out distortions. The Bi ions 
exist in the charge disproportionate state which is chem- 
ically interpreted as 2Bi''^+ ^^Bi'^^-f Bi^+. It is therefore 
tempting to connect the mechanism of superconductivity 
with the nature of these insulators. 

Unfortunately, experimental estimates of the electron- 
phonon coupling strength do not lead to the firm 
conclusion on the origin of superconductivity in the 
bismuthates. Large isotope effect with a=Q.A has 
been reported for BKBOa. Other measurementda give 



a=0.21±0.03 for BKBO and a=0.22±0.03 for BPBO. 
Using their analysis, the authors of Ref. H concluded 
that "phononic" effects in these materials are only in- 
dicative of dressed electronic excitations. From studying 
the imaginary part of optical conductivity in BKBO the 
authors of Ref. ^ gave the value of A ~0.2. Electronic spe- 
cific heat measurementsQ of Bao.6Ko.4Bi03 have yielded 
N*{0) ~0.32 states/ [spin xeVx cell] giving a mass en- 
hancement factor N*{0)/N^'''"^{0) -1.4 (7V,^""''(0) -0.23 
states/ [spin xeVx cell] for x=0.4). A complicated sit- 
uation exists with the transport measurements. The 
temperature dependent resistivities for superconducting 
BKBO and BPBO have ranged from metallic to semi- 
conducting and two-channel model of the conductivity 
in the bismuthates was discussedEl. While good grain- 
boundary-free thin films and single crystals of BKBO 
doped well away from the CDW instability seem to ex- 
hibit metallic behaviour, the values of the resistivity itself 
are (like in the cuprates) unusually high and are of the 
order a few hundred /ifixcm at room temperature This 
could point out that an additional (to standard electron- 
phonon) scattering mechanism is presented. 

The most direct evidence on the importance of 
electron-phonon interactions in superconductivity of 
the bismuthELtea. has been given by the tunnelling 
measurementgj^EJ. Although not identical for differ- 
ent junctions, the deduced Eliashberg spectral functions 
a'^F{uj) bear a close resemblance with the phonon density 
of states determined by inelastic neutron scatteringllil. 
The estimated values of A vary from 0.7 to 1.2 which 
seem to be sufficient to explain high critical temperatures 
within the standard mechanism. 
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The electron-phonon coupling in doped BaBiOa has 
been investigated theoreticaUy by several methods. The 
authors of Ref . |l2| study this problem using tight -binding 
(TB) fit to the energy bands which are obtained from 
density-functional calculations based on local-density 
approximatior£3 (LDA) . The computed Eliashberg spec- 
tral function a^F(uj) has been found to display promi- 
nent features in the frequency range corresponding to 
the oxygen stretching modes and the value of A=1.09 
has been reported. Crude calculations based on rigid- 
mufhn-tin approximation (RMTA) raJso give large A ~3 
indicating a strong-coupling regimctil. Two estimates of 
A using total-rejiergy frozen-phonon method appeared in 
the literaturetSllS. Note that, in contrast to the TB and 
RMTA methods, the frozen-phonon calculations treat 
screening of the potential due to lattice distortions self- 
consistently. The value of the electron-phonon coupling 
strength equal to 0.3 for the breathing mode has been. 
foundE£l and the rough estimate of A '^0.5 was obtainedEEl 
using 12 q=0 phonons for ordered cubic Bao.sKo.sBiOs. 

There was a partial success in predicting structural 
phase diagraHiriipr the parent compound BaBiOa within 
the LDAllaO"t3. The experimental structure mainly 
consists of combined tilting and breathing distortions of 
the oxygen octahedra corresponding-.to the instable R- 
point phonons of the cubic phaseE^O. While rotational 
instability was found in all calculations, the frozen-in 
breathing laode was not described by pseudopotential 
calculationtLajid two linear-mufRn-tin-orbital (LMTO) 
calculationsEjO give the value for the breathing distor- 
tion about 30% off the experimental one. Less rigor- 
ous potential-induced-breathing (PIB) model obtained 
both instabilitiesE3 with similar accuracy. It is not clear 
whether these discrepancies are due to sensitivity to com- 
putational details or due to the local density approxima- 
tion itself. p 

A great amount of worko has been done to understand 
the properties of the barium bismuthates on the basis of 
the negative U extended Hubbard model originally intro- 
duced in Ref. The valence configuration of semicon- 
ducting BaBiOs can be viewed as Ba2Bi'^+Bi^+06 which 
represents a lattice of electron pairs centred at every sec- 
ond Bi site (Bi'^+). The sites occupied with Bi^+ ions 
are intexoreted as those with no electrons. Rice and co- 
workerfj have proposed that such local pairs are sta- 
bilised by polarising the O octahedra and the efi^ective 
on-site U becomes negative due tOjthe large electron- 
phonon coupling. Recently, VarmacJ has pointed out 
that negative U can also be of the electronic origin due to 
the skipping of the valence "4-f " by the Bi ion. The latter 
can provide a possible explanation for a well-separated 
optical and transport energy gap in the bismuthateso. 
The mean-field phase diagram of the negative U Hub- 
bard model exhibits several stable phases involving a 
CDW semiconductor, and a singlet superconductor. This 
is in qualitative agreement with the experimental phase 
diagramcZl. 

The question on the origin of negative U is of great 



interest since it may provide an insight on the supercon- 
ductivity mechanism in the bismuthates. Recent calcu- 
lations using constrained density-functional theory have 
been carried out to obtain the-Coulomb interaction pa- 
rameters for the Bi 6s orbitalf£3. No indication for neg- 
ative U of the electronic origin was reported. 

In the present work we try to address several prob- 
lems seen from the above introduction by means of the 
state-of-the-art density functional LDA calculations. As 
a first problem, we study lattice dynamics of the super- 
conducting cubic Bao.6Ko.4Bi03 using recently developed 
linear-response approach implemented within the LMTO 
mcthodcj. The effect of doping is modelled by the vir- 
tual crystal and virtual mass approximations. On the 
basis of this calculation, we estimate electron-phonon 
coupling in this compound. The linear-response method 
used by us is advantageous in contrast to the frozen- 
phonon approach since it allows the treatment of pertur- 
bations with arbitrary wave vectors q. We have demon- 
strated its accuracy by calculating lattice-dynamical, su- 
perconductin^and transport properties for a large vari- 
ety of metalsEJ, and we believe that our calculated value 
of A will be a realistic estimate for the electron-phonon 
coupling strength in-this high-Tc superconductor. Also, 
a recent publicationcll deals with the application of the 
linear-response method to study the electron-phonon in- 
teraction in another high-Tc superconductor CaCu02. 

The second problem, we focus in our work, is study- 
ing the effects of anharmonicity in the electron-phonon 
coupling. It is widely accepted that certain phonon 
modes are strongly anharmonic in the high-Tc mate- 
rials. Frozen-phonon calculations produce double-well 
potentials for buckling motions of oxygen atoms_j)ei:- 
pendicular to the CuO planes in nearly all HTSCEiTLa, 
chain-buckling distortions are found to be anharmonicEj 
in YBa2Cu307, X-point tiltings of the octahedra along 
(110) directions are instableEJ in La2Cu04, and R-point 
instabilities corresponding to breathing and tilting exist 
in the doped barium bismuthatest£Hl3£3. The influence 
of anharmonicity to high-temperature supepiO|tLductiv- 
ity has been addressed in several publicationsOtZl, espe- 
cially because of the small isotope efi^ect found for HTSC 
cuprates. Triple-degenerate pure rotational mode at the 
R point of cubic ordered Bao.sKo.sBiOa was predicted 
to exhibit a double-well potential and some estiniates of 
the anharmonic contributions to A have been givenlij. We 
extend this analysis by solving numerically Schrodinger's 
equation for the anharmonic potential well found from 
frozen-phonon calculations. The anharmonic A is then 
computed along the lines proposed in Refs. by es- 

timating the electron-phonon matrix elements from the 
energy bands computed for different tilting distautions. 
We conclude, in accord with the previous findingsll3, that 
this contribution, while not decisive, is not negligible for 
the total value of A. 

The third purpose of our work is to study the prop- 
erties of the undoped parent compound BaBiOa. We 
try to answer the question whether the LDA gives an 
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adequate description of the ground state properties for 
this charge-density-wave insulator. Since there jsras-some 
inconsistency reported in previous calculationsEj~t3, we 
want to rule out possible sensitivity of the final results to 
the internal parameters used in our band structure cal- 
culations with the full-potential LMTO methocO. We 
carefully choose our LMTO basis set, number of k points, 
plane wave energy cutoff and other parameters by exam- 
ining the convergency of the total energy and the cal- 
culated properties with respect to them. Based on the 
well converged data, we come to the conclusion that the 
breathing distortions are seriously underestimated (ide- 
ally, absent) in the LDA, and, therefore, the insulated 
state is not correctly described. This strongly resembles 
the situation with the antiferromagnetic ground state of 
the cuprates-superconductors which is also not described 
by the LDAc2l. We perform a number of model calcula- 
tions in the spirit of the LDA-I-U methodcil in order to 
clarify this problem. 

The rest of the paper is organised as follows: In Sec. II, 
our linear-response calculations of the lattice dynamics 
and the electron-phonon interactions in doped BaBiOa 
are described. Sec. Ill considers anharmonicity correc- 
tions to A for the tilting motions of the oxygen octahedra. 
Sec. IV reports our LDA and model calculations of the 
ground state properties for pure BaBiOs. In Sec. IV we 
give our conclusions. 

II. HARMONIC PHONONS AND A 

This section presents our results on the lattice dy- 
namics and the electron-phonon interaction for the cu- 
bic perovskite superconductor Bao.6Ko.4Bi03. We also 
summarise the main features of the calculated electronic 
structure and discuss our predicted equilibrium lattice 
configuration. The band structure calculations are per- 
formed pwi^h the highly precise full-potential LMTO 
methodo. The details of the calculations are the fol- 
lowing: The effect of potassium doping is taking into 
account using virtual crystal approximation (VGA) by 
considering a fractional nuclei charge Z=5^.6 at the Ba 
site. Numerous supercell investigational3'c3 of the dop- 
ing influence on the calculated energy bands justify the 
applicability of the VGA. A multiple, three-K LMTO ba- 
sis set with the tail energies equal to -0.1,-0.8, and -2 Ry. 
is employed for representing valence wave functions. The 
valence states include 6s and 6p orbitals of Bi, 2p orbitals 
of O, and 6s orbitals of Ba. Such semicore states as 5d 
orbitals of Bi, 2s orbitals of O, and 5p orbitals of Ba are 
treated as bands and arc included in the main valence 
panel using two k LMTO basis with 2=-0.1,-0.8 Ry. 
The main panel also includes unoccupied 5d orbitals of 
Ba with the 2k basis and 4/ orbitals of Ba with the Ik 
basis (k^=-0.1 Ry). Deeper lying 5s states of Ba are re- 
solved in a separated energy panel. All other states are 
treated as core levels. The mufHn-tin sphere radii were 
taken to be: S'Ba=3.25 a.u., S'Bi=2.25 a.u., and 50=1.80 
a.u. All calculations are performed at the experimen- 
tal lattice constant a=8.10 a.u. The Barth-Hedin-likc 



exchange-correlation formula after Ref. 43 is used. The 
valence bands are treated scalar relativistically and the 
core levels - fully relativistically. A number of k points 
for the Brillouin zonis-iBZ) integration using an improved 
tetrahedron mcthodca is taken to be 20 per-gth BZ. The 
charge density and the potential in the interstitial region 
are expanded in plane waves with the cutoff correspond- 
ing to the (28,28,28) fast-Fourier-transform (FFT) grid 
in the real space (approximately 10000 plane waves). 

We first summarise the main features of the calculated 
electronic structure in doped BaBiOs. The occupied part 
of the bands (see Fig. |^) mainly consists of Bi(6s)-0(2p) 
hybridised band psmplex. This is in accord with the pre- 
vious calculationEj. For the cubic perovskite phase, there 
is only one band crossing the Fermi level, which is an an- 
tibonding Bi-0 sp{a) band. Similar situation is found 
in the cuprate superconductors where Cu-0 dp{a) an- 
tibonding bands dominate at the Fermi energy Ep. A 
simple tight-binding model involving Bi(6s), 0{2p) or- 
bitals, and two-centre nearest-neighbour sp{a) interac- 
tion can be used tn_understand the principal features of 
these energy bandsc3. It was early notedc3 that for the 
case of half-filling (undoped cubic BaBiOa) this model 
has a perfectly nested Fermi surface for the wave vector 
q corresponding to the R-point. Therefore, it is tempt- 
ing to interpret the appearance of breathing distortions 
as commensurate Peierls instability and cubic perovskite 
Bao.6Ko.4Bi03 as a doped Peierls insulatorEj. To under- 
stand whether nesting can bring any effect in static sus- 
ceptibility, we have analysed q~dependence of the func- 
tion 
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FIG. 1. Calculated LMTO energy bands for cubic 
Bao.eKo. 46103. The potassium doping is taken into ac- 
count using virtual crystal approximation. 

for the realistic energy bands i^kj (relative Ep) us- 
ing experimental structures. We conclude that the nest- 
ing is far from perfect in the case of half-filling and 
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dimerisation of the oxygen octahedra can hardly be con- 
nected with it. A reahstic TB model should also in- 
clude Bi(6p) orbitals and— their nearest-neighbour inter- 
action with 0{2p) statesota. Upon potassium doping, 
the bands hardly change except for a slight lowering of 
Ep away from half-filling. For Bao.6Ko.4Bi03 the Fermi 
surface represents a rounded cube centred at the T point 
as shown in Fig. H. Analysis of the band structure fac- 
tor given by Eq. ffl) as a function of q shows featureless 
behaviour and any effect of the nesting enhancement on 
the electron-phonon interaction is not expected for this 
band dispersion. 




phase. According to the neutron diffraction data@, the 
average structure is cubic although presence of a weak 
long-range superstructure characterised by the octaJ; 
dra rotations at the angles about 3° was also foundol 
Recent XAFS measurementsEJ report on the locally dis- 
ordered rotations. From their analysis, the authors 
of Ref. ^ conclude that the rotations can either be 
along (1,1,1)- or (1,1,0) axe. Previous frozen-phonon 
calculationgli3 performed for the ordered Bao.sKg.sBiOa 
investigate (1,0,0) component of the tilting mode which 
is found to be unstable with the total energy minimum 
corresponding to the angle 7°. 

Our own total-energy calculations also confirm the ex- 
istence of tilting distortions. Fig. ^ shows that the total 
energy exhibits a double well behaviour as a function of 
the rotation angle. We choose (1,1,0) axe for the tilting 
as is the case in the undoped compound. The unit cell 
in the calculation is doubled according to the R point of 
the cubic phase. The total energy minimum is found at 
the angle equal to 5°. The energy gain compared to the 
cubic phase is only 10 meV/(lxcell) which indicates that 
at the temperatures of the order Tc the rotations can be 
dynamic. The double well behaviour at such small en- 
ergy scale unambiguously points out on the importance of 
evaluating anharmonicity contribution in total clcctron- 
phonon coupling. This problem will be discussed in the 
following section. 



FIG. 2. Calculated Fermi surface for cubic Bao.6Ko.4 
BiOs. using the LMTO method. The effect of potassium 
doping is taken into account within the virtual crystal ap- 
proximation. The centre of the cube corresponds to the F 
point of the Brillouin zone. 

We second discuss our results for the calculated equi- 
librium lattice configuration in Bao.6Ko.4Bi03. The 
thcoretical-to-experimental volume ratio V/Vexp is 
found to be 1.01, and the calculated bulk raadulus is 
equal to 1.25 Mbar. Both neiitron diffractiorO and x- 
ray-absorption-fine-structurecZI (XAFS) measurements 
show that frozen-in breathing distortions are absent in 
the superconducting phase. We have performed frozen- 
phonon calculations for the doubled cell corresponding to 
the i?-point and for several breathing distortions. The 
total energy minimum shows that the undistorted cubic 
phase is stable in accord with these experiments. The 
curvature is well fitted with standard parabola, which 
shows that the breathing mode is harmonic in the super- 
conducting phase. 

We further investigate tilting of the octahedra. Ex- 
perimentallvr-, for the undoped compound the octahe- 
dra rotatedEll at the angle ~11.2° along (1,1,0) axe. 
More complicated situation exists in the superconducting 
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FIG. 3. Frozen-phonon calculation of the total en- 
ergy (meV/lxcell) as a function of the tilting angle in 
Bao.6Ko.4Bi03. The levels e„ are the solutions of the 
Schrodinger equation for the anharmonic oscillator with 
the double-potential well shown on the figure. The transi- 
tions a;„ = e„ — eo involving different phonon excited states 
are illustrated by arrows. 

We now report our main results on the calculated lat- 
tice dynamical properties of Bao.eKo. 46103. The den- 
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sity functional linear-response approached implemented 
on the basis of the full-potential LMTO methodu is used 
in this calculation. The dynamical matrix is computed 
at 20 irreducible q points corresponding to the (6,6,6) 
reciprocal lattice grid of the cubic BZ. The effect of the 
potassium substitution on the phonon spectrum is taken 
into account by virtual mass approximation. The LMTO 
basis set and other technique details have been described 
above. One more comment should be said on evaluating 
BZ integrals in the linear-response calculation. Here, 
one can essentially improve the accwacy of the integra- 
tion by using a multigrid techniquccJ. A (6,6,6) grid (20 
irreducible k-points) is used for finding linear-response 
functions while the effects of the energy bands and the 
Fermi surface are taken into account using a (30,30,30) 
grid (816 irreducible k-points). 

The calculated phonon spectrum along major symme- 
try directions of the cubic BZ is plotted in Fig. ^ . Solid 
circles denote the calculated points and the lines result 
from interpolation between the circles. Several features 
can be seen from these phonon dispersions. Three high- 
frequency optical branches around the lo ^^17 THz are 
well separated from the other modes distributed in the 



frequency range from to 10 THz. The high-frequency 
modes mainly consist of the oxygen bond-stretching vi- 
brations. The longitudinal branch at the point R cor- 
responds to the breathing mode which in our calcula- 
tion has a frequency 15.7 THz. From the analysis of 
our polarisation vectors, we conclude that oxygen bond- 
bending vibrations dominate in the frequency interval 
between 6 and 10 THz. The octahedra tilting modes 
are at the low-frequency interval. They exhibit signifi- 
cant softening near the q-point i?=(l,l,l)7r/a. Due to 
symmetry, one can talk about pure tilting at the line 
between q-point A'/=(l,l,0)7r/a and the R point. Ex- 
actly at the M point nearest octahedra tilt in-phase and 
they tilt out-of-phase at the R-point. A nearly-zero- 
frequency triple-degenerate mode exists at q=(l,l,l)7r/a 
which corresponds to the pure rotational T2u phonon. In 
fact, for T—Q this mode should have a slightly imagi- 
nary frequency for the cubic structure according to our 
frozen-phonon analysis illustrated in Fig. ^ . But, due 
to numerical inaccuracies the linear-response calculation 
gives very small positive a;=0.5 THz. No significant soft- 
ening of the tilting modes near the point M is predicted 
by our calculation. 
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FIG. 4. Calculated phonon spectrum of Bao.sKo.dBiOs using density-functional linear -response method. The 
potassium doping is taken into account using virtual crystal and virtual mass approximation. The calculated points 
are shown by symbols. The lines result from interpolation between the points. Horizontal bars indicate the measuredc3 
phonon frequencies. Numbers for every phonon mode indicate the calculated electron-phonon couphng constants Aqi^. 
(Only the values larger than 0.01 are emphasised.) On top of the figure shown are (i) the values of A summed over all 
branches for given q, (ii) the values of q-dependent transport constant Atr. The calculated phonon density of states 
F{uj) is shown on the right. 
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The phonon dispersion curves along TX, TM, and TR 
symmetry directions for Bao.6Ko.4Bi03 have been . very 
recently measured by inelastic neutron scatteringcSlEa. 
Horizontal lines in Fig. ^ indicate measured phonon fre- 
quencies at the symmetry points T, X, M, and Re as we 
were able to deduce them from Fig. 1 of Ref. ^ The 
existence of soft rotational modes near R can be directly 
seen from the measured phonon dispersions. The authors 
of Ref. 49 have reported that their samples still have a 



week long-range superstructure characterised by the tilt- 
ing of the octahedra. An extremely low frequency ~0.9 
THz of these modes was measured. This is in agreement 
with our calculations. 

Two other comments should be said on the comparison 
between our theory and the experiment. One comment 
concerns frequency interval from to 10 THz. Here, our 
calculation is seen to reproduce measured phonon disper- 
sions with the accuracy of the order 10%. In particular, 
the lowest mode in F has ujcaic— 3.79 THz which can be 
compared with Uexp "^3.5 THz. This mode mainly in- 
volves Ba(K) and Bi vibrations. The next mode in F is 
the oxygen out-of-phase mode. Here, u)caic=i-88 THz 
and the measured frequency is less than 6 THz. The 
so-called ferroelectric mode has a frequency 5.59 THz in 
our calculation which is close to ojexp found near 6 THz. 
This mode is bond-bending longitudinal and it has the 
strongest polar character. Usually it exhibits large split- 
ting from the TO mode at F in cubic perovskites. The 
presence of free charge carriers screens Coulomb interac- 
tion at long distances, and therefore, the LO-TO split- 
ting is absent in our calculation. The dispersion of the 
ferroelectric mode as a function of q is also seen to be 
correctly reproduced. 

The second comment concerns our comparison for the 
high-frequency interval, where the results of the calcu- 
lations are found to be less accurate and the overall dis- 
crepancy consists about 20%. The highest mode at F is 
Bi-0 bond-stretching mode. Here, we report the value of 
ll> equal to 17.91 THz and the oj value only slightly larger 
than 15 THz was found experimentally. The authors of 
Ref. ^ discuss an anomalous dispersion for the longi- 
tudinal optical branch of the one-dimensional breathing 
mode along FX with its pronounced frequency renormal- 
isation. Our calculation, on the other hand, gives much 
less dispersive optical branches along this direction as can 
be seen from Fig. ^ It is not clear whether this result 
is due to inhomogeneity of the potassium distribution or 
other imperfectness of the samples used in the experiment 
or due to drawbacks in our calculation connected with 
the virtual crystal approximation. In fact, it is clearly 
seen that all our high-frequency branches are overesti- 
mated by ~20% in contrast to the experimental ones (ex- 
cept eventually the breathing vibrations near the point 
R). This result also follows fpam the comparison of our 
calculated and the measuredU phonon density oL state 
^"(0;) [see FigJ^(a)]. It was found experimentallyliil that 
the oxygen bond-stretching modes exhibit softening with 
the substitution of Ba by K. These modes are located at 



the energies ~70meV (or 17 THz) in undoped BaBiOa. 
Therefore it is tempting to connect possible source for the 
discrepancies with our poor treatment of doping. The 
authors of Ref. |l^ discuss doping induced appearance of 
localised holes on the oxygen 2p orbitals which screen 
the charge on the oxygen anions. This charge reduction 
will lower the energy of these modes. If the localised hole 
picture is correct, the VGA will not capture this since it 
removes electrons from the conduction band by(uniform 
distributing the holes between 0{2p) and Bi(6s) orbitals. 

We now report our results for the calculated electron- 
phonon interaction. Based on our screened potentials 
which are induced by nuclei displacements and are found 
self-consistently, we evaluate matrix elements of the 
electron-phonon interaction, q?'1 -,, -. The standard 
expressioncil for the electron-phonon matrix elements 
reads as 



.9k+qj'kj 
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where Q^^' are the orthonormalised polarisation vectors 
associated with the mode qiy, Mr are the atomic masses, 
R runs over basis atoms in the unit cell and a runs over 
directions x,y,z; 5^^V denotes self-consistent change in 
the potential associated with the q-wave displacements 
of atoms R along a axe. In practical calculations we have 
also added so-called incomplete basis set corrections to 
the matrix elements (^) according to the method devel- 
oped in Ref. |o| 

The coupling strength Aq^ for the electrons with the 
phonon of wave vector q and branch j/ is given by the 
following integral 
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where Ns{Q) is the density of states at Ef=0 per cell and 
per one spin. Indexes j and j' numerate the bands (not 
spins) and spin degenerate case is assumed throughout 
the paper. The total coupling constant A results by sum- 
ming Aqj/ over ly and by averaging over BZ. Two delta 
functions in (^ impose integration over the space curve 
resulting from the crossing of two Fermi surfaces sepa- 
rated by q. For this integral we have used as many as 
816 k-points per irreducible BZ. 

The calculated values of X^i, at the symmetry direc- 
tions of the BZ are indicated in Fig. ^ along with the 
calculated phonon dispersions. (We emphasise only the 
values larger than 0.01.) On top of the figure shown 
are the values of Aq which are summed over all branches 
for given q. It is seen that the electron-phonon cou- 
pling is large for the high-frequency bond-stretching lon- 
gitudinal branch. Tliis result is expectable from band 
structure calculationsc3 since bond-stretching and espe- 
cially breathing vibrations produce modest changes in 
the energy bands near Ep. Coupling is strongly en- 
hanced near the points M and R, where it reaches the 
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values ~0.3. Here, the mode corresponds to two- or 
three-dimensional breathing. The value of Ah=0.3 for 
the breathing mode at R is-in accord with the previous 
frozen-phonon calculationEj. The authors of Ref. ^ on 
the other hand, give much lower value for equal to 
0.04. For other bond-stretching vibrations, we find Xqi, 
of the order 0.1. 

From Fig. ^ we conclude that the elect ron-phonon 
coupling is not small for the bond-bending oxygen 
modes. It is seen that A is enhanced for the wave vec- 
tors near R and also along TX direction. In the latter 
case, the vibrations correspond to the ferroelectric mode 
and the value of Aqj^ as high as 0.43 is found for the q- 
point (l/3,0,0)7r/a. Strongly anharmonic tilting modes, 
on the other hand, do not exhibit noticeable electron- 
phonon coupling in the linear order with respect to the 
displacements. Exactly at the point R, these modes have 
electron-phonon matrix element equal to zero by sym- 
metry, and, therefore, small values of to do not bring any 
effect on enhancing the coupling. We refer to the follow- 
ing section on our evaluated anharmonicity corrections. 




0.6 
0.5 



modes ^ 




30 40 50 60 70 

Frequency, meV 



80 



FIG. 5. Results for doped BaBiOpri (a) Comparison 
between calculated and experimentalcil(symbols) phonon 
density of states. (b) Calculated Eliashberg spec- 
tral function ret^F{Lu) and the results of the tunnelling 
measurementsB (symbols), (c) Calculated transport spec- 
tral function a trF{uj). 



On the basis of our evaluated phonon dispersion and 
q-dependent electron-phonon interaction we calculate 
the Eliashberg spectral function a^F{Lu). This is plot- 
ted in Fig. ^(b) by full lines. There, we also show 
by symbols two a'^F{uj) -which were deduced from the 
tunnelling measurementsO. Comparing the experiment 
with our calculations, it is first seen that the intensities 
of high-energy peaks are approximately the same which 
means that we reproduce the coupling for these phonons 
sufficiently accurate. It is also seen both from the the- 
ory and the experiment that while in the phonon density 
of states [Fig. ^(a)] mainly the TO phonons contribute 
to the high-energy structure, for the a^F(ijj) these are 
the LO phonon modes. Our tendency to overestimate 
the phonon frequencies at high energies is again clearly 
distinguishable . 

The most prominent feature seen from our calculated 
a'^F{uj) is the absence of any structure for the low- 
frequency interval below 40 meV. This strongly con- 
tradicts with the experimental a^F^uj) which exhibits 
two intensive peaks in this region centred at 15 and 30 
meV. According to our analysis of partial F{u;), the ori- 
gin of the first peak could be due to low energy Ba(K) 
and Bi vibrations together with the tilting modes, and 
the second peak can result from the bond-bending oxy- 
gen modes. It is not clear however why the calculation 
seriously underestimates the electron-phonon coupling 
for these phonons, while it correctly describes the cou- 
pling for the bond-^tiretching modes. From the band 
structure argumentsEa one can expect that only bond- 
stretching modes will have large interaction with elec- 
trons. The bond-bending modes cannot produce any 
significant changes in the bands near Ep since sp{a) in- 
teraction for these kind of distortions is not changed in 
linear order. The same is true for the tilting modes. Since 
there is no partial weight of the Ba(K) orbitals at Ep we 
also do not expect strong electron-phonon coupling for 
the low-frequency interval. 

Another possible explanation for the observed peaks 
is due to the contributions connected with anharmonic- 
ity. For the anharmonic phonons, not only one-phonon 
virtual transitions must be seen in the a^F(tj), but also 
higher order virtual states. The detailed discussion on 
this subject will be presented in the next section, here 
we only give the value 0.04 as our final answer for the 
anharmonic contribution to A resulting from the tilting 
motions. While not negligible, this value alone again does 
not explain the intensity of the experimental a'^F{uj) at 
low-energies. 

Our calculated total value of A resulting from the har- 
monic phonons and linear electron-phonon coupling is 
found to be 0.29. This is too small to account for the su- 
perconductivity at 30K in the compound BoosKoABiO^. 
In fact, only the high-energy phonons contribute to our 
coupling. As a result, our estimated value of wjog high 
as 550K is found. Using McMillan's expressionE3 with 
/i*=0 we however find the critical temperature with our 
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set of parameters as low as 4.5K. One can try to esti- 
mate the error in our A value due to the overestimation 
of the frequencies for the bond-stretching modes. Using 
the expressionEj: A ~ NP /Mq'^, where NP is the elec- 
tronic prefactor and MuP' is an average force constant, 
one sees that A should increase with lowering uP' . Our 
20% error in u) results in 30% error in <D^, and this can 
lead to the actual A values which are 30% higher than 
we calculate. However, A ^^0.4 is also not sufficient to 
explain the value of Tc- 

Our calculated elect ron-phonon contribution to the 
transport properties is a final subject of this section. The 
quantities responsible for the electronic transport a*€ eas- 
ily deduced from the linear-response calculationa23. By 
inserting the electron velocity factor (wkj — Wk+qj' f to 
the expression (|^), we calculate transport constant \tr- 
Its q-dependence is shown at the top of Fig. |4| along with 
the q-dependence of the electron-phonon A. We see that 
both functions exhibit very similar behaviour in the BZ. 
The transport spectral function alj.F{uj) is the central 
quantity needed for evaluating temperature-dependent 
electrical and thermal resistivity as lojsst-order variational 
solutions of the Boltzmann equatiorO. Our calculated 
a^^F{Lo) is shown in Fig. ||(c). It is seen that this 
function behaves closely to the supierconducting a^F(w) 
which is usually the case in metalst3. The total average 
Xtr is found to be 0.32. Based on these data, we evaluate 
electron-phonon limited electrical resistivity p to be 14.3 
jiD, X cm at T=273K. This is at least one order of magni- 
tude lower than the values of p reported in the literatureo. 
It is clear that the source for this discrepancy is the same 
as in our describing superconducting properties. It is un- 
likely that the strong electronic correlations are presented 
in these materials because the parent compounds are dia- 
magnetic (not antiferromagnetic) insulators. Therefore, 
it is unlikely that other (than electron-phonon) scatter- 
ing mechanisms take place, as spin fluctuations, for ex- 
ample, in HTSC cuprates. Taking into account anhar- 
monic phonons, polarons, or bipolarons may be decisive 
for describing the superconductivity and transport phe- 
nomena here. Our basic conclusion is that the conven- 
tional ideas on the electron-phonon mechanism are not 
operative in the hismuthates. 

III. TILTING AND ANHARMONIC A. 

Towards further understanding of superconductivity 
in the bismuthates, we try to evaluate anharmonicity 
corrections in the electron-phonon interaction. As we 
have mentioned in the introduction, several experiments 
and numerous frozen-phonon calculations point on pos- 
sible imttortance of these effects in the theory of HTSC 
cupratest^^EJ. In the bismuthates, the first candidate 
to study anharmonicity is the tilting mode: Our linear- 



response calculations give nearly zero-frequency vibra- 
tions for this pipde at the i?-point of cubic BZ. Our own 
(and previousll3) frozen-phonon calculations predict here 
a double-potential-well behaviour with very shallow en- 
ergy minimum at 5° as illustrated in Fig. |^. Despite the 
reported-average structure forj-the superconducting phase 
is cubicEll, some experimentsca discuss the existence of 
long-range superstructure characterised by the xotations 
of the octahedra. Recent XAFS measurementsEJ report 
on the presence of locally disordered rotations. 

Unfortunately, though formulated0, the problem of 
the influence of anharmonicity to superconductivity is 
not tractable numerically in a full volume. Our sim- 
plified treatment is. based on the expression introduced 
by Hui and Allenc3 which generalises zero-temperature 
electron-phonon coupling to the anharmonic case by in- 
cluding matrix elements over all phonon excited states. 
The phonon states |n) and their energies e„ are obtained 
by solving the Schrodinger equation for an oscillatory 
mode characterised by q and v. (We will not label the 
states |n),e„ with [c^v) for simplicity). For harmonic po- 
tential wells e„ — eg is just nuj^i, (in atomic units), where 
cjqjy is the phonon frequency. This leads only to the one- 
phonon virtual states {n^\) which arc involved in the 
matrix elements of the electron-phonon interaction. For 
anharmonic potential wells the spectrum is generally dif- 
ferent from the set of equidistant levels. One example is 
a double well of the tilting mode which is shown in Fig. 
^ The real spectrum e„ obtained as the solution of the 
Schrodinger equation is plotted in Fig. ^ by horizontal 
lines. Therefore, our first purpose is to examine what 
effect in our A would bring the proper treatment of all 
virtual phonon states. The second problem is connected 
with the modification of the electron-phonon matrix ele- 
ments due to higher-order terms in the expansion of the 
change in the ground state potential with respect to the 
atomic displacements. Since we wish to examine these 
effects only for the tilting mode at the point R of the 
cubic BZ, we make an essential approximation by as- 
suming that the tilting mode is not coupled to the other 
modes of either this wave vector or other wave vectors 
which is generally not true when anharmonic terms are 
included into the lattice dynamical problem. We assume 
that the polarisation vectors for this mode are known and 
are given by exactly the out-of-phase rotations of the 
nearest octahedra along the (1,1,0) axe. Any processes 
of phonon-phonon interactions will be neglected in this 
treatment. We also neglect by the finite-temperature ef- 
fects using the arguments given in Ref. ^ In the double- 
well problem one expects that the modifications due to 
frequency and electron-phonoH|pHriatrix-elements renor- 
malisation are quite dramaticr^rl and bring the largest 
effect in the values of A. 



We start from a general, zero-temperature expression for the electron-phonon coupling in anharmonic case, which 
can be written as foUowa^a 



8 



1 



2 



where w„ = e„ — cq. For the moment we will not assume as in Ref. ^ that uin are small at the electron energy scale. 
(The latter reduces the integral with the Fermi step functions /kj — /k'j' to the integral with the delta functions as 

given by Eq. ^.) We introduce a generalised matrix element 0^^^}^^^ for the virtual transition to the n's phonon state: 

G'K'k. = (^l^3k','k,|0). (5) 

The electron-phonon matrix element A^k'j'kj in (|^) involves the transitions between the states |kj) and |k'j') near 
the Fermi surface 

A^k'.'k, = (k'j'|Al^|kj-), (6) 

where AV is the total change in the ground state potential induced by the lattice distortion associated with the mode 
CIV. We have especially included A in the notation Af/k'j'kj since AF is not a derivative of the potential with respect 
to the displacements but it is the difference between the self-consistent potential V{y, {tR + AtR}) for the distorted 
crystal and the potential T^(r, {tR}) for the undistorted crystal. The atomic positions at the equilibrium are given by 
the vectors tR = t + R, where t denote the translations and R are the basis vectors. The displacements associated 
with the mode (\v are described by the vector field AtR. Therefore AV is proportional to AtR and so does A^k'j'kj- 
By introducing complex (infinitesimal) polarisation vectors 5Q_r of the modeE3 qi/, the displacements in any atomic 
cell t can be found using the formula: 

AtR„ = JQRc^e^'i* + c.c. (7) 

where a runs over directions x, y, z, and c.c. stands for the complex conjugate. (The quantities Ay, AtR, and (JQ/j 
should, in principle, be labelled with (\v but we omit this for simplicity.) 

The phonon states \n) are the functions of the displacements A^r^ or ^Qrq,. In order to compute the matrix 
element (n|A(7k'j'kj|0} over the phonon states we should expand AV^ with respect to the displacements. Keeping the 
terms up to second order, this expansion reads as 

Ay = ^<5QR„^e^'i* 



Ra t 



RR'qq' tt' 

Here, 5V/6t-£ia is associated with the first-order derivative of the potential when a single nucleus centred at t + R 
experiences an infinitesimal displacement along a-th direction, and / 5t-^a5t'^, ^, stands for the second-order 
derivative. [Notation iRo, is shorthand for (t + R)a.\ Obviously, both these response functions have no dependence on 
the mode (\v. If V{t) has a periodicity of the original lattice, the change 5V/5t-R,a is a function of general type. One 
expects that SV/Stna is only not zero in the vicinity of the displaced atom and it goes to zero when r departs from 
the site t + R. However, because of the translational invariance of the original crystal, considering the response at 
the point r due to the movement of atom in t + R must be equivalent to considering the response at the point r — t 
due to the movement of atom at R (when t=0). Therefore we can write that 5V{Y)/5tYia — 5V{t — t)/6Ra- We now 
introduce the lattice sum 

sLv = Y.^'-'^ (9) 

which represent a variation of the potential per unit displacement induced by to the movements of atoms R along 
Qf-th axe by infinitesimal amount St-^io, proportional to exp{iqt) in every cell t. It is easy to prove that the expression 
(^ translates like a Bloch wave with wave vector q in the original lattice, i.e. 6~^V{r + R) = e^'^^6~^V{r). (We 
will sometimes omit indexes Ra) Notation S~^V refers to the travelling wave of vector +q while complex conjugated 
quantity S^V would refer to the travelling wave of vector — q. 
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One can analogously define lattice sums associated with the second-order changes of the potential. These enter the 
second and third contributions in Consider, for example, the lattice sum associated with the second contribution: 



,5(2)1/ 



This expression translates like a Bloch wave of vector 2q because 



(10) 



E 



E 



<5(2)y(r) 



„iq(t+t') 



tt 



<5(i - t")Ra5(t' - t")n'. 
:q(t+t')j!!!n£L. 



(11) 



Analogously, the lattice sum associated with the third contribution in ^ can be denoted as 5~^S~V. It represents a 
travelling wave of wave vector 0^ i.e. it is periodical at original lattice. 

Using the notations and (|To|), the change in the potential AV given by the formula now has the form 



Rq 



RR'aa' 



(12) 



RR'i 



It is clear that when this expansion is used in the matrix element (j^), certain selection rule will occur for the wave 
vectors k' and k. Namely, the matrix element (k'j'|5+X^|kj) is equal to zero unless k' = k + q, (k' j'\5^S^V\]<.j)~0 
unless k' = k + 2q, and (k' j'\6+6-V\kj}^0 unless k' = k. 

Let us now introduce the electron-phonon matrix element associated with the first-order change in the potential 



•-^k'i'kj — Ok'k+q 



k + qj" 



(13) 



The electron-phonon matrix elements associated with the second-order changes in the potential have two forms 
according to the second and third contributions in ( p^ 



GL"i^k; = ^W.(k + 2qj' 



J2 {n\SQRo.SQR,c.'\0) X S+J+,^,V 



RR'i 



(14) 



^[n]{2'} _ 1 . / , ., 



kj 



(15) 



\ RR'qq' 

Then, the expression (^ for Xq^, splits into three contributions 

A = A^iJ + A^2} , ;^{2'} 

associated with one electron-phonon matrix element from the first order, Eq. (p^, and two matrix elements from the 
second-order, Eqs. (Q, and (|l5|), i.e. 



(16) 



a{i> 



A{2} 



A{2'} 



Ns{Q) 



E '^if^j ^ /k+qj')^(^kj - £^k+qj' + >^n) 



■-^k+qj'ki 



E E^-'^'^J ~ ./'k+2qj')'^(-^kj - -E^k+2qj' + W„) 



kjj' n 



^N{2} 
'-^k+2q/kj 



kjj' n 



G 



n]{2'} 
kj'kj 



(17) 
(18) 
(19) 
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The double sums over k and k' appeared in (m are now reduced to the single sums over k according to the selection 
rules in the matrix elements (|l3|), (^^, and (|l5|). Taking into account matrix elements coming from the complex 
conjugated quantities 6~V, 6~S~V, and S^S^V gives contributions with wave vectors k — q, k — 2q, and k. However 
making substitutions k — q ^ k and k — 2q ^ k in the k-space integrals of (||) results in extra factor of 2 which 
appears in (|l^), (p^), and ( [l9| ) compared to (^). 

We now discuss the derived expressions. First, it is instructive to see how Eq.|l^ for X^J goes to the standard 
formula (|^) . For the case of harmonic oscillator we can use the properties of Hermite polynomials and prove that the 
matrix element over the phonon states appeared in (|l3|) is reduced to 

(n|<5QR„|0)=<5„i-^^==, (20) 

where Qrq denote ortho nornialised polarisation vectors of the mode qy in contrast to JQrq which are related to 
the actual nuclei displacementa£j. The matrix element ( ^ ) is only not zero for the transitions including one-phonon 
virtual state. Substituting (20) into (13) gives exactly the formula (^. By placing ( p^ ) to ([l^) we see that only n=l 



term in the sum over n survives, and the energy difference ui is equal to ujqi^. The standard formula is then 
recoveredEil by replacing (/kj — f-k.ij')5{E\^j — £'k+qj' + ^JJc^u) to Wf^y5{E\^j)5{E\^^qj') which is valid for small uj^i,. 

Second, it should be noted that the number of matrix elements necessary to reach the convergency in the sum over 
n is actually not large. From the numerical estimates of the oscillator strengths /„ for the transitions from the phonon 
ground state |0) to excited states the authors of Ref. ^ concluded that the f-sum rule (sum over all oscillator 
strengths gives unity) is closely satisfied by taking into account only siugle term /i. Even for infinite rectangular well 
which represents an extreme case for anharmonicity, n=5 is sufficientE3. Our own numerical experiments with the 
double-potential well of the form shown in Fig. |^ confirms this conclusion both for the dipole matrix elements and 
for quadrupole ones. This implies that the phonon excitation energies w„ = e„ — eo appeared in the expressions (p^), 
(19|) are still too small at the electronic energy scale and, therefore, all integrals over k are reduced to the 



m, and 



integrals over the space curve resulting from the crossing of two Fermi surfaces separated by q. It then follows that 
the expression for is always equal to zero unless there are electronic inter band transitions with zero momentum 
transfer at the phonon energies. The latter standardly is not assumed. We therefore left with the following formulae 
for and Aq^-'^: 



A^l^ = ]^E'^(^'^^-)'^(^w)E|Gwk,r/^n, (21) 

^^'^ ^im^ c5(£;k,)'5(i^k+2ci,-0 E r (22) 

Another comment concerns the expressions (|l^) and ( p2[ ) for Aq^'^ which involves electronic intraband transitions 
with momentum transfer equal to 2q. If wave vector q is such zone-boundary vector that 2q = G, where G is one of 
the reciprocal lattice vectors, then it is obvious that Aq^- — according to the arguments given above. T he situ ation 



here is analogous to the discussion of contributions to A from the T optical phonons (see, for example Ref. ^,^6|). It is 
well known that exactly q=0-phonons do not couple to the electrons since the energy cannot be conserved. Therefore, 
small but finite vectors q LUqiy/vp, where vp is the Fermi velocity, must be considered. 



We now discuss implications of the theory given above 
for the tilting mode in Bao.eKo. 46103. Exactly for the q- 
point i?=(l,l,l)7r/a the contribution to A connected with 
the first-order changes in the potential, X^J , vanishes 
due to symmetry of the matrix elements. Also the con- 
tribution connected with the second-order changes, Aq,y , 
is zero since vector 2q is equivalent to the T point. The 

contribution Aqi/ is always zero since there is only one 
band crossing Ep and there is no interband transitions 
near Ep. In the higher orders, there are contributions 
to A connected with the odd-order changes in the poten- 
tial and with the even-order changes. In the odd orders, 
changes in the potential would represent travelling waves 



of vectors q, 3q, 5q, etc., and the contributions to A would 
be equal to zero due to symmetry of the matrix elements. 
In the even orders, changes in the potential would repre- 
sent travelling waves of vectors 0, 2q, 4q, etc., which are 
all equivalent to the T point when q=(l,l,l)7r/a, and cor- 
responding contributions to A also vanish. Therefore, we 
conclude that pure rotational mode at R does not couple 
to the electrons in any order. 

In order to obtain an estimate of the effect, we must, 
in principle, step out from the point R. Unfortunately, 
for general q we cannot perform frozen-phonon calcula- 
tions of anharmonic coefficients, and the problem looses 
its numerical tractability. Another approach is to dis- 
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cuss quantities resulting from the integration of the ex- 
pressions (H^), and (|2^ ) over q . In this way, we should 
assume that the matrix elements Gk+qj^k^ and G{^2q/kj 
have no q-dependence (which is generally not true but 
the order of magnitude of the effect will be certainly 
captured) and can be approximated by the values at 
the point R. Therefore, the contribution resulting from 
the integration of disappears due to symmetry of 
^k+<i/'kj ^"^^ integration of Aq'^^ over q gives: 



|r^[«]{2} 



(23) 



(Here we have assumed that k + 2q is equivalent to k.) 
The obtained expression is an average of squared defor- 
mation potential which is induced by the tilting of the 
octahcdra in the second order with respect to the dis- 
placemcnts.._Lt has the same meaning as discussed in the 
literatureH'EB magnitude of the change of the phonon 
self-energy due to transitions across the gap in the su- 
perconducting state (so called superconducting A^). 

The value of Xl ' can be readily evaluated using the 
frozen-phonon method. Let us measureE3 the tilting 
distortions by tilting angle t. Quadrupole matrix ele- 
ments between the phonon states are then given by 
(n|i^|0) and the electron-phonon matrix elements are 
(kJlS'^^^V/ 5t'^\kj) where S'^^^V/St'^ is the second-order 
derivative of the potential with respect to the tilting an- 
gle taken at t=0. Thejmatrix elements {kj\S^^^V/St'^\kj) 
can be approximatectj by the second-order derivatives 
of the one-electron energies i?kj relative to Ep , They are 
accessible from band structure calculations. Using this, 
the formula (Esl) for A^^-*^ can be rewritten as follows 



A^^>=25:^(i?k,)^^^^^x 



dt2 



E 



\{n\e\Q)\^ 



(24) 



where the quadrupole matrix elements and the energy 
differences uon are obtained by solving the Schrodinger 
equation for the double potential well shown in Fig. 
Only even-order matrix elements {n = 2k\t^\Q) are al- 
lowed by symmetry of our distortion, and therefore, the 
transition between nearly degenerate first and second, 
states shown in Fig. |^ (i.e. wi 0) does not contributeLj 
to Xi^\ 

In order to evaluate ([24|), we first calculate energy 
bands E\^j\t\ relative to Ep[t] for a set of distortions t, 
and then integrate 



F{t)^2Y,5{E^M)^{E^A 

kj 



(25) 



over the unperturbed Fermi surface. Due to symme- 
try, the energy bands i?kj [i] have only even terms in 
the Taylor's expansion over t around t=0. The function 
F{t) — F[Q) is fitted to the polynomial starting with at"'. 
Then, our \]^^ is 



2n 



^2n 



(26) 



Using frozen-phonon calculation, we evaluate a=2.991 
Ry/rad^. Solving the Schrodinger equation for the 
double well, we also evaluate the sum over n in ( |2^ ) 
equal to 0.0044 rad-^/Ry. This gives the total estimated 
Ai —0.013. Since the tilting mode is triple degenerate, 
the final result A'f^^ = iX]^^ ^ 0.04 is obtained. It can 
be compared with_the value 0.11 previously reported by 
Kunc and ZeyhertS who used a somewhat analogous ap- 
proach. Unfortunately no details were presented in Ref. 
|l6| and, therefore, it is difficult to determine the main 
source of our discrepancies. 

Despite several approximations were made in the above 
derivation, we think that the values of the order 0.04 
(to 0.11) can be realistic for anharmonic contribution 
to A from the tilting modes involving large ionic excur- 
sions. This constitutes about 20% of the A value 0.29 
found by the linear-response method. Unfortunately, 
our total harmonic plus anharmonic Atot=0.33 is still 
too small to account for the superconductivity at 30K 
in Bao.6Ko.4Bi03. 

In order to have large effect in A due to anharmonic- 
ity, one has to analyse three contributions in Eq. |2j. 
The first one is connected with the changes in the energy 
bands. The second contribution goes from the matrix 
elements |(n|t2|0)p which can be large if large atomic 
displacements are possible due to the flatness of the po- 
tential well of the oscillator. This is exactly the case for 
tilting distortions. The third contribution can go from 
the energy denorairiator in case values ti;„=e„ — eq become 
sufficiently smallEll. 

The fact that our anharmonic A'^^'^ is small is due 
to small changes in the one-electron energies associated 
with tilting. This is so because the tilting distortions do 
not change the distance between Bi and O atoms, thus 
keeping sp{a) nearest-neighbour interaction nearly con- 
stant. Possible large contribution from the energy bands 
would be expected for breathing distortions. Unfortu- 
nately we have not found strong anharmonicity for the 
breathing potential well in doped BaBiOa. It is how- 
ever known that in the undoped compound there exist 
frozen breathing distortions which correspond to a deep- 
double-well situation. It could be possible that our LDA 
description of breathing distortions is not completely cor- 
rect. The following section is devoted to this question. 



IV. BREATHING AND LDA 



So far we have discussed the compound Bao.6Ko.4Bi03 
which is a metal in its normal state. Our calculations 
show that such ground state properties as equilibrium 
structure and lattice dynamics in the adiabatic approx- 
imation are reproduced by the density functional LDA 
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tice. In ad- 
The nature 



method reasonably well. This suggests that our descrip- 
tion of the elect ron-phonon interactions based on the 
LDA energy bands is realistic. We now turn to the dis- 
cussion of the undoped parent compound BaBiOa. It is 
known for many years that pure BaBiOs is a charge- 
density-wave semiconductor in a sense that two bismuth 
atoms exist in the charge disproportionate state Bi''^^. 
This leads to the modulated Bi-0 distances (breath- 
ing distortions) in the cubic perovskite 
dition, there are strong tilting distortions 
of the disproportionate state is still not very well under- 
stood fiUHjd several explanations involving Fermi surface 
nestin^J, a real sjiace pairing based on strong electron- 
phonon interactioncJ, and the existence-pf negative elec- 
tronic U due to Bi^'^ valence skippingE3 have been sug- 
gested in the past. It is also unclear whether there is any 
connection between the origin of the semiconducting be- 
haviour and the doping induced superconductivity a the 
border of metal-insulator transition. 

If intra atomic correlations of Bi 6s electrons are strong 
and responsible for the appearance of negative U, one 
expects that the LDA theory will fail to describe the 
charge-disproportionate state. This exactly happens in 
HTSC cuprates where due to large positive U of Cu 3d 
electrons, an antiferromagnetic ground state of undoped 
cuprates was not predicted by the LDAcJ. If, on the 
other haxid, negative U is due to strong electron-lattice 
couplingcil, the LDA should quantitatively explain the 
observed instabilities. (It is likely that the nesting idea 
is not very convincing as we discussed in Sec. II of this 
paper. The breathing instability, for example, has not 
been found in La2Cu04 where the same logic is valid.) 

Several total energy LDA calculations of the struc- 
tural phase— .diagram for pure BaBiOa exist in the 
literaturetj'EIrt3 and the results |-seem to be inconsis- 
tent. Pseudopotential calculationgl3 predicted rotational 
but not a breathing instability. In contrast to that, 
pure breathing distortion b ~ 0.06A with large energy 
lowering AiJf, ~ -50 meV/(lxcell) has been founclE£l 
by the linear-augmented-pclane-wave (LAPW) calcula- 
tion. A third caJculationEJ based on the full-potential 
LMTO methodEJ gives here nearly zero AEf, ^ —0.7 
meV/(lxcell) with the equilibrium value-i) ~ O.OSA for 
pure breathing. Using the same methodu, the authors 
of Ref. m calculate AEb = -20 meV/(lxcell) and b = 
0.07A. The calculations of the total energy Etot as a func- 
tion of combined tilting plus breathing (tb) distortions 
have also been reported. The first such calculationtj 
predicts AEtb =-14 meV/(lxcell) relative to rotational 
energy mirppum, t = 8.5° and b ~ 0.05 5 A. The second 
calculationEl gives AEtt =-40 meV/(l xcell), t = 9.6° 
and 6 = 0.11 A. Note that the experimental values found 
for r=150K are: texp = 11-2° and &e^p = 0.085A. AJfga 
rigorous potential-induced-breathing model predictsc3'Ej 
these properties with similar accuracy. In addition, the 
reportedllj frequency Wb=46 meV of th e br eathing mode 
is too low comparing to the measuredOca values ^ 70 
meV. 



These large discrepancies in predicting structural prop- 
erties of BaBiOs seem to contradict with the accuracy of 
LDA which standardly is of the order a few per cent. This 
could already signalise that the ground state of this semi- 
conductor is not completely captured within mean-field 
LDA solution. Also, another possible reason for the ob- 
tained discrepancies could be due to unusual sensitivity 
of the calculated properties to the computational details. 

We have performed our own studies of the structural 
phase diagram for pure BaBiOa basejd on the highly- 
precise full-potential LMTO methodLa. We have indeed 
found that there is a sensitivity of the final results to 
the details of the calculations, and we will discuss this 
due course. However, we have also found that serious er- 
rors are introduced by the LDA. Our general set-up for 
the calculations is the same as was used for the doped 
compound Bao.6Ko.4Bi03. It has been described at the 
beginning of Sec. II of this paper. 

As a first step we check an equilibrium cell volume. 
The value of V/Vexp — 0.998 is found and the calculated 
bulk modulus is equal to 1.29 Mbar. These values es- 
sentially depend on the treatment of the semicore states. 
Especially, treating deep lying 5d Bi states in the main 
valence panel, i.e. allowing their full hybridisation with 
the valence states, is found to be crucial for the appear- 
ance of the total energy minimum itself. Such sensitivity 
is due to the fact that the ground state potential V^(r) 
is highly not spherical. In the Bi-0 directions, the po- 
tential exhibits a pronounced minimum (~ 0.5 Ry down 
its average value at the Bi MT~sphere) , and the average 
kinetic energy of Bi 5d orbitals in the interstitial region 
is nearly zero. As a consequence, the character of Bi 
5d electrons appears at O atoms. In the Bi-Ba direc- 
tions the potential exhibits large maximum (^^ 0.5 Ry up 
its average value at the Bi MT-sphere), and the kinetic 
energy here is large negative. Kinetic energy variation 
about 1 Ry leads to necessity to use multiple k LMTO 
basis even for semicore states. 

We now report our results for the calculated tilting in- 
stability. The unit cell is doubled in this calculation ac- 
cording to the -R-point of the cubic phase. The rotations 
of the octahedra arexierformed along the direction (1,1,0) 
as in the experimentEJ. The calculated Etot versus tilting 
angle is shown in Fig. ^. A pronounced minimum corre- 
sponding to the angle t = 13° is found by our calculation, 
and large energy lowering AEt = —200 meV/(l xcell) is 
predicted. Note that for the doped compound, the calcu- 
lated energy lowering becomes only -10 meV/(l xcell) as 
can be seen from Fig. ||. The tilting angle calculated by 
us agrees with the expcrimcnicj within the same accu- 
racy as previously reportcdij'ta. However, our calculated 
AEt disagrees with AEt^-24: meV/(l xoeU) of Ref.Ej 
and with AEt=-60 meV/(l xcell) of Ref.H. To check 
the sensitivity of our values to the treatment of semicore 
states, we have performed the calculations by placing Bi 
5d orbitals into a separate energy panel. This indeed has 
an effect in the equilibrium tilting angle equal now 9° 
and A£'t=— 50 mcV/(l xcell). Four times change in the 
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latter value looks very unusual and again shows the im- 
portance of proper handling with the semicore. We have 
also tried to increase the number of k-points in the inte- 
gration over the BZ but this hardly has any effect to the 
final results. 
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FIG. 6. Frozen-phonon calculation of the total energy 
(meV/lxcell) as a function of the tilting angle for undoped 
BaBiOa. 

Wc now discuss our studies for the breathing distor- 
tions. The dimerisations of the octahedra are performed 
for a set of tilting angles from to 15° in order to search 
for a global total-^energy minimum at the bt plane. The 
crystalline structure consists of two formula units and 
has now a monoclinic symmetry. Our calculations do 
not predict pure breathing distortions {t = 0°) leaving 
cubic structure stable against this perturbation. How- 
ever, for t — 0° our total energy displays very flat highly 
anharmonic potential well assuming closeness to the in- 
stability. When tilting appears, the potential well starts 
to exhibit a double minimum indicating weak but non- 
zero breathing distortion. We have found that global 
total-energy minimum occurs when t = 13° The calcu- 
lated Etot versus b for this angle is shown in Fig. |^. 
A very shallow minimum can be seen from this figure 
located at 5 ^ 0.04A. The corresponding energy low- 
ering AEi, is only equal to -7 meV/(lxcell) relative to 
the structure which is purely tilted. These repalts are 
in contradictions with the experimental findingscJ which 
give bexp = 0.085A. Moreover such low AEi, gives us the 
breathing phonon frequency nearl y equal to zero in con- 
trast to the measured valuestilcS of the order 70 meV. 
This indicates that real energy lowering is much larger 
than we evaluate. 

Furthermore, the LDA one-electron spectrum calcu- 
lated at our minimum still corresponds to metallic ground 
state. We have found that the energy gap Eg opens at 



larger values of 6 0.07A. At b = b^xp = 0.085^ the 
calculated minimal gap is indirect and occurs between 
the points X and L (see also Ref. 



15). Its value is ap- 



proximately equal to 0.1 eV. The minimal direct gap is 
found at the point L and is about 1 cV. These valuta 
are lower than the corresponding experimental valuesE£l 
most likely due to systematic underestimation of gaps 
by the LDA.. The measured transport activation gap is 
0.24 eV. However, it does not show up in photoconductiv- 
ity, optical absorption or photo- acoustic measurements. 
Sometimes this is interpretedED as bosonic bound state 
of two electrons due to negative electronic U . There- 
fore it is unclear whether transport activation gap can 
be related to the minimal indirect gap of our calculation. 
The gap seen from the optical measurements is about 
2 eV which is substantially larger than our 1 eV direct 
gap. We have however not performed the calculations 
of the optical properties, therefore an exact comparison 
between the theory and experiment is not currently pos- 
sible. 
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FIG. 7. Frozen-phonon calculation of the total energy 
(meV/lxcell) as a function of the breathing distortion (A) 
for undoped BaBiOa. The upper curve is obtained using 
68 k-points for the integration over |BZ of the monoclinic 
lattice and the lower curve corresponds to 476 k-points. 

Several test calculations have been made to check our 
results. This mostly concerns the breathing instability, 
the energy bands were found to be insensitive to the com- 
putational details. First, we place Bi bd orbitals into a 
separate energy panel. This has an interesting effect that 
the breathing become more pronounced, the calculated 
value of 6 = 0.065A becomes much closer to the experi- 
ment, and the energy lowering is now -17 meV/(lxcell) 
relative to the purely tilted structure. This result, how- 
ever, should be considered as artificial, since separate 
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treatment of Bi 5d states does not allow us to reproduce 
the equilibrium cell volume. We second tried to analyse 
the convergency with respect to the multiple n LMTO 
basis. We have used up to 5 k basis functions for repre- 
senting valence wave functions but this practically does 
not affect the final results. Adding higher lying Ba 6s and 
Bi 6d orbitals has no effect either. This indicates that the 
basis set described above (see Sec. II) is sufficiently com- 
plete. We third try to investigate the relativistic effects 
due to heavy Bi atoms. Inclusion of spin-orbit coupling 
along with the scalar relativistic terms was found to have 
a negligible effect at the calculated equilibrium structure. 
Fourth, we increased the number of k-points to 476 per 
|;th BZ of the monoclinic lattice. This has an effect of 
lowering total energy for metallic states (i.e. when b < 
0.07A) and practically does not change the values of Etot 
for semiconducting energy bands. As a result, total en- 
ergy as a function of b becomes even more shallow (see 
Fig. ^ and predicted breathing at the equilibrium is 
extremely small. 

Based on our findings we conclude that the breathing 
distortions are seriously underestimated (ideally absent) 
in the LDA, the predicted ground state is metallic and, 
therefore, the charge-density-wave instability is not cor- 
rectly described. This situation is analogous to that with 
HTSC cuprates where antiferroHjagnetic ground state 
was also not found within LDAe3. (It therefore seems 
that LDA may, generally, have problems with both spin- 
density waves and charge-density waves.) We also think 
that due to either improper handling with Bi 5d orbitals 
or approximate treatment of the full-potential ternjia.ia. 
the LMTO method of Ref. ^ previous calculationaljy 
did not converge to the true LDA ground stale. 

Following to the Hohenberg-Kohn theoremE^ we know 
that more proper treatment of exchange-correlation ef- 
fects beyond LDA should in principle reproduce semicon- 
ducting ground state of BaBiOs and the correct values 
for breathing distortions. At the same time, the energy 
gap will not necessarily be reproduced since it is not a 
ground state property of a single system. However, in 
such systems as BaBiOs the energy gap is directly re- 
lated to the charge disproportionation between two Bi 
atoms since the splitting between occupied and empty 
Bi 6s levels is proportional to the occupancies of these 
orbitals. The latter is related to the charge density dis- 
tribution which is a ground state property. It is therefore 
seems that until the correct theory will not reproduce the 
energy gap values, the correct breathing distortions will 
not be obtained. Speculating on this, (i) we do not see 
how the LDA can describe the correct ground state and, 
at the same time, strongly underestimate the gap value, 
(ii) we think the exact density functional theory would 
describe both the ground state and the energy gap in this 
system. 

It worth to mention in this context an example with 
antiferromagnetic oxides like NiO. It is known that LDA 
calculations underestimate both the magnetic moment 
and the energy gap value in this compound. It is also 



clear that the magnetic moment is given by the occu- 
pancies of 3dx2_y2 and 3dz2 orbitals of Ni, and the same 
factor defines the splitting between occupied and empty 
states which is directly related to the energy gap. A so- 
called LDA plus U mcthocH provides more proper treat- 
ment of the systerris-|With strong electronic correlations. 
It has been shownEJ that this method predicts correct 
ground state for many Mott-Hubbard insulators and, at 
the same time, gives more accurate values for the energy 
gaps comparing to the LDA. Also, similar improvements 
are obtained with— the use of self-interaction corrected 
density functionalE^. If we would interpret these results 
as using a better energy functional, then the situation 
with NiO should be completely analogous to BaBiOs 
with the exception that not the spin transfer but the 
charge transfer defines the properties of this system. 

These conclusions may bring an attention to the exis- 
tence of intra atomic correlations of Bi 6s electrons not 
captured by the LDA. It does not necessarily follows from 
our calculations that the Coulomb U is negative in these 
systems. What only follows from our calculations is that 
there exists a correction term to the LDA in which the pa- 
rameter responsible for the attraction of two electrons at 
the Bi sites is negative. To illustrate this, let us restrict 
ourselves for the moment by only Bi 6s orbitals. The cor- 
rection energy for the doubled unit cell Ba2Bi^+'Bi^~'06 
can be represented in the form (see Appendix) 

A£:eorr = ^ ^^'^ffi'^^ " = ^^efjq'^ (27) 

i=Bil,Bi2 

where Ui is the occupancy of the Bil=Bi''+'^ or 
Bi2=Bi*~' 6s state and n = (nsii +'^Bi2)/2. The charge 
disproportionation parameter q = {ubh — 'T._Bi2)/2. When 
AUeff < 0, the correction ( p7| ) to the LDA total en- 
ergy is negative for any non-zero q, and therefore the 
charge disproportionate state becomes favourable. This 
construction is purely heuristic and isJpuilt on the anal- 
ogy to the LDA-I-U density functionalcil. (A better name 
here could be "LDA minus U".) In order to avoid double 
counting effects we interpret AUeff as U — Ulda, where 
Ulda is a part of the on-site interaction taken into ac- 
count in the LDA. (In this way, when the LDA is ade- 
quate for describing the correlations, the correction term 
becomes automatically zero.) Therefore, from our calcu- 
lations follows that the LDA overestimates Coulomb U 
for Bi 6s electrons, and the difference AUeff is less than 
zero. Adding AEcorr to the LDA functional will obvi- 
ously result in obtaining the charge disproportionation. 

Unfortunately this simplified illustrational model has 
several drawbacks. First, the antibonding band crossing 
the Fermi level (see Fig. |l|) consists not only of the Bi 
6s electrons but also has a substantial character of the 
O 2p and Bi 6p electrons. Therefore, AUeff in Eq. ( [27| ) 
should be interpreted as the on-site interaction between 
the correspondingly constructed Wanier functions. Sec- 
ond, since these Wanier functions are long ranged, inter- 
site Coulomb interaction parameters must be introduced 
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into the expression (p7|). (In fact, by Fourier transform- 
ing the antibonding band at Ep we evaluate the range 
of Wanier states to be 4 lattice constants of the cubic 
phase.) Another possibihty is to consider a multiband 
model involving charge fluctuations and Coulomb inter- 
actions between Bi 6s, O 2p and Bi 6p states. 

The Coulomb interaction parameters can, in princi- 
ple, be I calculated using constrained density functional 
methodE3l£j. However, the number of parameters re- 
quired for our system seems to be much larger than just 
one number AJ7e/ / • This gives a lot of extra freedom and 
complicates their finding. Despite this difficulty— some 
progress has already been made in this directionnS. We 
also plan to investigate the parameters required for the 
extended multiband Hubbard model using constrained 
LDA calculations. This question will be addressed in the 
future publications. 

In the following, we take the simplified model ( p7| ) to il- 
lustrate which features will bring the inclusion of the cor- 
rection energy to the LDA functional. Using a variational 
principle to minimise our LDA-C/ expression for the to- 
tal energy Eld A + AEcorr leads to solving the single- 
particle equations with the potential Vlda + ^Vcorr: 
where AV^^^^ = AUeff{ni — n) = AUeffq is the con- 
tribution which has an orbital dependence. It is clear 
that the value 2AV^g„ determines the effective splitting 
between two non-equivalent Bi 6s levels. If our LDA cal- 
culation gives this splitting nearly equal to zero (unless 
we set the breathing distortion to the experimental one), 
then its real value should be of the order of optical gap 
experimentally determined as 2 eV. (Simple connection 
between the charge disproportionation and the gap value 
is clearly seen from here, and we again state, that we 
do not see the way how the LDA can obtain the correct 
breathing and at the same time seriously underestimate 

We perform the calculations iavolving AV^^^^ using 
our full-potential LMTO methocO. The calculations are 
analogous to the constrained LDA calculations described 
in Ref. |6^. Using projector operator technique we define 
the corrections to the Bi 6s diagonal matrix elements of 
the LMTO Hamiltonian. The projection is simply taken 
to the I = spherical harmonic inside the MT sphere of 
the Bi site and zero everywhere else. According to our 
model, the difference between two non-equivalent Bi 6s 
occupation numbers should be associated with the charge 
disproportionation 2q. If we now fix AUef / to some value, 
the self-consistent procedure will define this charge dis- 
proportionation, and the correction AEcorr to the LDA 
total energy can be estimated. In this way, we can find 
a new equilibrium structure and compare that with the 
experiment. Choosing such AUeff to reproduce the ex- 
perimental breathing distortion can give some insight on 
this value. 

We have performed a set of self-consistent calcula- 
tions of the total energy Elda + AEcorr as a function 
of breathing distortion for different values of AUeff- We 
have indeed found that for AUeff of the order minus 10 



eV, the experimental breathing distortion is fairly well 
reproduced and the correct breathing phonon frequency 
is obtained. Also, the energy gap becomes much closer 
to the experiment. It is interesting to note that the 
form (27) of the functional assumes broken symmetry, i.e. 
the existence of the charge disproportionation regardless 
the presence of breathing distortions. As a consequence, 
the double-potential well does not exhibit a smooth be- 
haviour at 5=0 and has a kink there. 

The value of AUeff =-10 eV required to describe the 
experimental structure seems to be unphysically large. 
This is mainly due to small self-consistent values of q of 
the order 0.1 electrons which are obtained in our calcula- 
tions. (Even smaller values ^ O.Ole are obtained for the 
charge transfer between two Bi MT-spheres because of 
the screening by the tails of O 2p orbitals). In fact, if we 
take a simple relationship 2qAUeff — 2AVcorr Eg — 2 
eV, the same AUeff =-10 eV is obtained. This again 
shows that our logic is valid. As we discussed above, a 
better way would be to interpret 2q as a charge transfer 
between two Wanier states centred in the sites of cu- 
bic lattice with two non-equivalent Bi atoms. (We as- 
sociate these Wanier states with the Bi atoms, but it 
is clear that they represent a mixture of a whole set of 
orbitals.) Indeed, our tight -binding calculations for a 
one-band model with long-range interactions (fitted to 
describe the dispersion of the band crossing the Fermi 
level) show that introducing a small splitting between 
the effective Bi levels gives the charge transfer q of the 
order 1 e. (The chemical interpretation of the charge 
disproportionation 2Bi*+ — »Bi'^+-|-Bi^+ is thus valid for 
the Wanier functions and is, of course, wrong for the 
real charge densities.) Unfortunately, implementing the 
Wanier representation in our LDA-C/ functional is not 
straightforward and we postpone this for the future work. 
However, a rough estimate of AUeff operative between 
the Wanier states can be found by setting g to 1 in the 
relationship 2qAUeff ~ 2 eV. This will result in more re- 
alistic value of AC/e//=-l eV. Such values are expectable 
from the errors introduced by the LDA, for example, due 
to self-interaction effects. 

This discussion shows that the proper treatment of 
the parameters of the model (27) is required. It is also 
clear that inclusion of the intersite Coulomb interactions 
is necessary to obtain a quantitative description of the 
ground state in BaBiOa. The above calculations cannot 
be considered as estimates of U. Therefore we would await 
of making the conclusion that negative on-site U of Bi 6s 
orbitals is of the electronic origin until we will not eval- 
uate the errors introducing by the LDA. In the above 
discussion we were trying to argue that the insulating 
behaviour and breathing distortions in this system are di- 
rectly related to the charge disproportionation, and the 
LDA must fail in describing this ground state already 
because it is not able to reproduce the energy gap. On 
the other hand, the LDA+U like techniques were seen to 
be very perspective for performing accurate calculations, 
and we will try to address this subject in the future work. 
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V. CONCLUSION. 

In the present paper we have reported our density func- 
tional LDA studies of the compound Bai_a;K2:Bi03. For 
its superconducting phase (x—OA) we have performed 
the calculations of full wave-vector-dependent lattice 
dynamical properties and the elect ron-phonon interac- 
tions. These calculations were based on recently de- 
veloped linear-response approach implemented in the 
framework of the full-potential LMTO method. The fol- 
lowing conclusions are made on the basis of our stud- 
ies: The calculated phonon dispersion curves along ma- 
jor symmetry directions agree reasonably well with the 
results of the neutron-scattering experiments. Some dis- 
crepancies have been found to occur in reproducing the 
bond-stretching oxygen modes. They were attributed 
to the virtual crystal approximation used for the treat- 
ment of potassium doping. It was found that the bond- 
stretching modes have large coupling to the electrons. 
Especially the breathing modes at the R and M points 
of the cubic Brillouin zone have large coupling equal to 
0.3. Relatively strong coupling is also predicted for the 
bond-bending oxygen modes. However, our calculated 
averaged value of A was found to be 0.29 which is too 
small to account for the high-temperature superconduc- 
tivity in the doped barium bismuthates. These results 
were supported by the calculated transport properties 
such, e.g., as electron-phonon limited electrical resistiv- 
ity. 

Our own and previous frozen-phonon calculations 
predicted highly anharmonic double-potential well be- 
haviour for the tilting of oxygen octahedra correspond- 
ing to the point R. We have performed detailed stud- 
ies on the influence of anharmonicity corrections to 
the electron-phonon coupling. By neglecting the pro- 
cesses of phonon-phonon interactions and utilising zero- 
temperature treatment of Hui and AUcnEEl, we have 
worked out the formula for anharmonic contribution to A 
up to the second-order with respect to the displacements. 
It was found that anharmonic A can be not small if large 
ionic excursions take place as in the case of the tilting 
vibrations. Using the full-potential LMTO method and 
frozen-phonon approach, we have estimated that con- 
tribution equal to 0.04. Our total A is thus 0.33. We 
concluded that while not negligible, anharmonicity cor- 
rections due to tiling modes do not help to explain the 
superconductivity in the bismuthates, and therefore, its 
origin still remains open and intriguing problem. 

As a final issue, we have done the calculations of the 
structural phase diagram for the undoped parent com- 
pound BaBiOs. A low-temperature experimental struc- 
ture consists of combined breathing plus tilting distor- 
tions. In agreement with the previous conclusions, we 
have found that tilting distortions are reasonably well 
reproduced by the LDA calculation. However, we have 
also found that the breathing distortions are seriously 
underestimated (if not absent) within the LDA. This 
contradicts with the previous findings most likely due 



to improper handling with the semicore states. The un- 
derestimation of breathing leads to the predicted ground 
state which is metallic while it should be insulating of 
the charge-density- wave type. This situation closely re- 
sembles antiferromagnetism of HTSC cuprates which was 
also not predicted by the LDA. Using a simple correc- 
tional scheme in the spirit of the LDA+U method we 
have tried to argue that the problem of breathing is 
most likely due to underestimation of energy gap by the 
LDA. While we cannot make any definite conclusions on 
whether negative U in the bismuthates is of the electronic 
origin, the failure of the LDA indicates that the elec- 
tron correlations of the Bi 6s electrons are not properly 
treated. Since, LDA calculations contain all electron- 
lattice coupling effects, we conclude that in order to re- 
cover an insulating state, an existence of the correction 
to the LDA with some negative AUef / of the electronic 
origin must be assumed. We believe that LDA±J7 like ap- 
proaches will help in further understanding of the physics 
in these systems from density functional point of view. 
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Appendix. 

The form of the correction can be understood using 
the following arguments: Suppose the LDA total energy 
Eld A gives non-degenerate ground state described by 
the density po as is the case of predicting BaBiOa without 
breathing distortions. Let us expand El]ja{p) around 
its minimum. Due to extremal property, this expansion 
starts from the second-order variations: 
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If the predicted ground state is wrong, then the first 
candidate which can take responsibility for this is the 
second-order derivative. The value of 6'^'^^ Elda/ 6p5p at 
p = Po is positively defined since Elda{p) has a min- 
imum here. On the other hand, the true density func- 
tional at Po would indeed have a local maximum or, more 
generally, a saddle point, since the true ground state is 
double degenerate corresponding to either breathing-in 
and breathing-out distortions or vice versa. Therefore, 
6^'^'' Etruel ^P^P should be negatively defined at p = po. 
(We assume that the first -order variation S'^^^ Etme / 6 p is 
equal to zero at p = po due to symmetry of the saddle 
point.) It is thus tempting to construct a density func- 
tional with the corrected second-order variation, i.e. 
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Assuming that kinetic energies of both Etme and 
Elda functionals are the same, the difference between 
the second-order derivatives is described by effective 
Coulomb interactions. This leads to the form (|2^) of 
the correction. 
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